J/ip suppression in the threshold model at RHIC and LHC energy 
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In the QGP based threshold model 0, 0], m addition to the normal nuclear absorption, J ftp's 
are subjected to an 'anomalous' suppression such that above a threshold density rij/^, all the J ftp's 
are melted. In the threshold model we have analysed the recent PHENIX data on the centrality 
dependence of J ftp suppression in Au+Au collisions at RHIC. Feedback from decay of the state \ 
is accounted for. J ftp's are anomalously suppressed above a threshold density, nji^ = 3.57 ± 0.17 
fm~ 2 . Threshold density for anomalous suppression of the state \ i s uncertain to a large extent, 
n x — 0.32 ± 0.32 fm~ 2 . The fraction F of the state % can n °t k> e determined unambiguously, 
depending on the nuclear absorption, it can vary from 20% to 40%. We have also predicted for the 
suppression in Pb+Pb collisions at LHC energy. In central Pb+Pb collisions, J ftp's are suppressed 
by a factor of 3-4. Suppression pattern is rather similar to that in Au+Au collisions, if not slighty 
less in central collisions. Using the PHENIX data on the participant number dependence of the 
Bjorken energy density, we have also estimated the QGP formation time. For critical temperature 
T c =192 MeV, estimated QGP formation time ranges between 0.07-0.09 fm/c. 

PACS numbers: PACS numbers: 25.75.-q, 25.75.Dw 
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I. INTRODUCTION 

In relativistic heavy ion collisions J ftp suppression has 
been recognized as an important tool to identify the pos- 
sible phase transition to quark-gluon plasma. Because 
of the large mass of the charm quarks, cc pairs are pro- 
duced on a short time scale. Their tight binding also 
makes them immune to final state interactions. Their 
evolution probes the state of matter in the early stage 
of the collisions. Matsui and Satz predicted that in 
presence of quark-gluon plasma (QGP), binding of a cc 
pair into a J ftp meson will be hindered, leadingto the so 
called J ftp suppression in heavy ion collisions [3j . Over 
the years, several groups have measured the J ftp yield 
in heavy ion collisions (for a review of the data prior to 
RHIC energy collisions, and the interpretations see Refs. 
[J, |5J]). In brief, experimental data do show suppression. 
However, this could be attributed to the conventional 
nuclear absorption, also present in pA collisions. 

In recent Au+Au collisions at RHIC, one observe 
a dramatic suppression of hadrons with high momen- 
tum, transverse to beam direction (high px suppression) 
0, 0, H, 0- This has been interpreted as an evidence 
for the creation of high density, color opaque medium 
of deconfined quarks and gluons [l(| • It is expected 
that high density, color opaque medium will leave its im- 
print on J ftp production. At RHIC energy, it has been 
argued that rather than suppression, charmoniums will 
be enhanced [ll|, [l2[ • Due to large initial energy, large 
number of cc pairs will be produced in initial hard scat- 
terings. Recombination of cc can occur enhancing the 
charmonium production. Recently, PHENIX collabora- 
tion have published their measurement of the central- 



ity dependence oiJ/tp suppression in Au+Au collisions 
at RHIC energy Data are taken at mid-rapidity 

< .35) and at forward rapidity (1.2 < y < 2.2). 
J ftp's are more suppressed at forward rapidity than at 
mid rapidity. In most central Au+Au collisions, while at 
mid rapidity, J ftp's are suppressed by a factor of ~ 3, 
suppression factor is ~ 6 at forward rapidity. PHENIX 
collaboration also measured the centrality dependence of 
J I tp su ppr ession in mid and forward rapidity Cu+Cu col- 
lisions 14]. Final analysis is not completed yet. Prelim- 
inary analysis indicate that, unlike in Au+Au collisions, 
J ftp suppression in Cu+Cu collisions, donot show large 
dependence on rapidity. Within the experimental errors, 
in most central Cu+Cu collisions, both at mid rapidity 
and forward rapidity. J ftp's are suppressed by a factor 
of ~ 2. Both the PHENIX data on J ftp production in 
Au+Au and in Cu+Cu collisions, are not consistent with 
models which predict Jf tp enhancement [HI [13] . It was 
also seen that various models, e.g . comover model [To] ], 
statistical coalescence model [16| or the kinetic model 
[i~7l . OH , fail to explain the (preliminary) PHENIX data 
on the nuclear modification factor for J ftp in Cu+Cu and 
in Au+Au collisions. The data are also not explained 
in the Glauber model of normal nuclear absorption p^ ]. 
The QCD based nuclear absorption model [2(| [2l[ could 
explain the PHENIX data on the centrality dependence 
of J ftp suppression in Cu+Cu collisions at RHIC but 
failed for Au+Au collisions [22} . Gunji et al [23j anal- 
ysed the PHENIX data in a hydro+J/^ model. The 
experimental J ftp suppression pattern in mid-rapidity is 
well explained by sequential melting of x, ip 1 and J ftp in 
an dynamically expanding fluid. Melting temperatures of 
J ftp and x are found to be Tj^ = 2.02T C , T x = 1.22T f 
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in agreement with lattice based model calculations [24 
They estimated the fraction of higher states as ~ 30% 
which is also in reasonable a gree ment with measurements 
in hadron-hadron collisions [25l . 1261 ] . 
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Recently, in the threshold model [HQ, we have anal- 
ysed the PHENIX data on J/tp suppression [H [H[. 
Blaizot et al [2, Q , proposed the threshold model to ex- 
plain the NA50 data [29| on anomalous J/tp suppres- 
sion in 158 AGeV Pb+Pb collisions at SPS energy. To 
mimic the onset of deconfining phase transition above a 
critical energy density and subsequent melting of J/ip's, 
J suppression was linked with the local transverse 
density. If the local transverse density at the point 
where J/tp is formed, exceeds a threshold value (nj^), 
J/tp y s are melted. Our analysis indicate that J/tp sup- 
pression in mid-rapidity Au+Au collisions are well ex- 
plained in the threshold model for threshold density 
n J/tp — 3.73 ± 0.29/m -2 [28j. Forward rapidity data 
also require similar threshold density but the threshold 
density has to be smeared greatly. Since large smearing 
of the threshold density is contrary to the spirit of the 
model, it was concluded that while J/tp suppression in 
Au+Au collisions at mid-rapidity are consistent with the 
threshold model, that in the forward rapidity are not. In 
[H], with the help of the PHENIX data on centrality de- 
pendence of Bjorken energy density, it was argued that 
the threshold density is related to the Bjorken energy 
density times the QGP formation time. Threshold den- 
sity nj/jp « 3.73/m~ 2 then indicate that QGP formation 
time is ~ 0.06-0.08 fm for critical temperature 192 MeV 
or less. 

The analysis in [2?], HH ignore two important aspect 
of J/tp production (i) feedback from higher states \ an d 
ip' and (ii) transverse expansion. Threshold model is a 
static model and transverse expansion can not be incor- 
porated in the model. However, the model can be ex- 
tended to include the feedback from the higher states. In 
the present paper, we have extended the threshold model 
to include the feedback from higher states and analysed 
the PHENIX data on mid-rapidity J/tp suppression in 
Au+Au collisions. We also use the threshold model to 
predict for the centrality dependence of J r /tp suppression 
in Pb+Pb collisions at LHC energy (/i=5500 GeV). 
The model predicts that the centrality dependence of 
J/tp suppression in Pb+Pb collisions is similar to that 
in Au+Au collisions at RHIC. The plan of the paper is 
as follows: in section II, we briefly describe the threshold 
model. PHENIX data on the centrality dependence of 
J / ip suppression in Au+Au collisions at mid-rapidity are 
analysed in section III. We have used the model to predict 
the centrality dependence of J/ip suppression in Pb+Pb 
collisions at LHC energy (section IV). In section V, we 
discuss the connection of threshold model with Bjorken 
energy density times the QGP formation time. Lastly, 
summary and conclusions are drawn in section VI. 



II. THE THRESHOLD MODEL FOR J/ip 
SUPPRESSION 

The details of the threshold model could be found in 
[HQ- It is assumed that fate of a J/tp depend on the 



local energy density, which is proportional to participant 
density. If the energy density or equivalently, the partici- 
pant density, exceeds a critical or threshold value, decon- 
fined matter is formed and all the J/tp's are completely 
destroyed (anomalous suppression). This anomalous sup- 
pression is in addition to the "conventional nuclear ab- 
sorption" . In conventional nuclear absorption, J/tp are 
destroyed due to inelastic J/ip- nucleon collisions. Trans- 
verse expansion of the system is neglected. It is implicitly 
assumed that J/ip's are suppressed before the transverse 
expansion sets in. 

In the threshold model, number of J/ip mesons, pro- 
duced in a AA collision, at impact parameter b can be 
written as, 

^(b) = J d 2 sT^(s)T^(b - s) 

x S anom (b, s) , (1) 
where T e ff(b) is the effective nuclear thickness, 

/oo poo 
dzp(b,z)exp{-a abs / dz/p(b, zl)), 
-oo J z 

(2) 

Cabs being the J/f/'-Nucleon absorption cross-section. 
Sanomfo, s) in Eq[T] is the anomalous suppression factor 
introduced by Blaizot et al. [1, 2]. Assuming that all the 
J/ip's get suppressed above a threshold density (rtj/^), 
the anomalous suppression can be written as, 

S an om{b,s) = @{nj/^ -n(b,s)) (3) 

It was observed [l[ that by smearing the threshold den- 
sity by a small amount, one can obtain better fit to the 
data, but at the expense of an additional parameter (A), 

S onom (b,s) = ~ [l - tanhA(n(b,s) - n J/v ,)] (4) 

Critical ingredient of the threshold model is the " local" 
transverse density. At impact parameter b and at the 
transverse position s, local transverse density it can be 
obtained as, 

n(b,s) = T A (s)[l-exp(-a NN T B (s-b))} 

+T B (b - s)[l - exp(-a NN T A {s))] (5) 

where un n is the nucleon- nucleon inelastic cross-section. 
At RHIC energy we use cat at =44 mb. As mentioned 
previously, the model give excellent description of the 
PHENIX data on participant number dependence of J / ip 
suppression with (i)a bs_~4.39mb and (ii) threshold den- 
sity nj/^ «3.73 fm- 2 [28|]. 

The model can be easily extended to include the feed- 
back from higher states x an d V^- Approximately 30% 
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of J ftp's are from decay of the state x and 10% from 
the state ip/ [25, 26] . Lattice model calculations indicate 
that dissociation temperature for \ aim are similar. 
T x w 1.2T C and Tf sa 1.1T C [2j[. J/tp's can survive 
much higher temperature Tj/^ ~ 2.1T C . We ignore the 
small difference between dissociation temperatures of \ 
and ip' and assume that fraction F of J ftp are from de- 
cay of the higher state \. Both the state J ftp and x can 



be anomalously suppressed. And like the state J ftp, the 
state x could also be suppressed due to inelastic collisions 
with the nucleons. x-nucleon absorption cross-section 
( a abs) nee d n °t be same as that for J/i/j-nucleon absorp- 
tion cross-section. If fraction F of the observed J / ip are 
from decay of the higher state x, the J ftp production at 
impact parameter b can be written as, 



j/i, 

AA 



(b) 



J/4, 

NN 



(1-F) 



d 2 sT 



eff (s)T° ff (b 



J 



s) x S_(b, s)+F I d 2 sT'f \ S )T"^ J (b - s) x S a x „ om (b, s) 



reff, 



(6) 



where the first term is as usual for the direct J ftp pro- 
duction. The 2nd term accounts for the feedback from 
the state x- T' e ^ is the effective nuclear thickness for 
the state x> 

dzp{b,z)exp(-a^ ba / dztp(b,z/)), 

-oo J z 

(7) 

and, S* nom is the anomalous suppression factor for the 
state X; 



S%nom( h , s ) = 9(n x -n(b,s)) 



(8) 



with threshold density n x . The function can again be 
smeared at the expense of an additional parameter A x , 



S* nom (b,s) = 1 [1 - tanhA x (n(b,s) - n x )] (9) 
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III. J/ip SUPPRESSION IN AU+AU 
COLLISIONS 



Threshold model as presented here has quite a few 
number of parameters, nji^, and n x , the threshold densi- 
ties for anomalous suppression of J ftp and X;the fraction 
F of the state Xj &abs and a* bs the J/-0-nucleon and 
X-nucleon absorption cross sections. Two more parame- 
ters, A and A x are added if smeared Theta functions are 
used for the anomalous suppression of J /ip and x- To re- 
duce the number of parameters, we make the simplifying 
assumptions that J/^-nucleon or x _n ucleon absorption 
cross section is proportional to square of the radius of 
the J ftp or the x meson, 



T X = 
abs 



J/4> 



&abs 



(10) 



Radius of different charmonium states has been calcu- 
lated, rj/^=0.25 fm and r x =0.36 fm [24[. The parameter 
a abs * s then expressed in terms of a a bs, cr*l = 2.07er a fc s . 



FIG. 1: Anomalous suppression as a function of transverse 
density, for threshold density 71^^=3.73 fm~ 2 and n x =0.37 
fm~ 2 . The solid and dashed line corresponds to F—0.1 and 
F=0.3 respectively. 



We also fix the smearing of the Theta functions to 
fixed values, A = 0.05tijm, X x = 0.05n x . Since dis- 
sociation temperature of the state x is approximately 
half of the dissociation temperature of the state J/ip 
[24j and in the threshold model, threshold densities 
are proportional to energy density, it is expected that 
the threshold density for the state x wn l be approxi- 
mately one-tenth of the threshold density for the state 
J/ip. In Fig[T] we have shown the anomalous suppres- 
sion S = (1 — F)Sanom + F 'S$ nom for a representative 
threshold density nj/^,=3.73 fm~ 2 and n x =0.37 fni~ 2 . 
The smoothened theta functions retain the basic feature 
of the threshold model, i.e. sudden onset of anomalous 
suppression of J/ip and x above threshold density ^3.73 
fm and ~ 0.37 fm~ 2 respectively. 
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FIG. 2: contours of constant^ 2 /-^ in a a bs — F plane. The 
threshold densities are is fixed at nj/^=Z.Q j'mT 2 and 
n x =0.36 fm~ 2 . One observe that equivalent X 1 /N is ob- 
tained by increasing F and decreasing a a bs or the vice-versa. 



With the simplifying assumptions, number of undeter- 
mined parameters in the threshold model is four, Jiju, 
n x , F and cr a bs- However, the parameters a a b s and F 
can not be determined simultaneously from a fit to the 
PHENIX data on the centrality dependence of J/ip in 
Au+Au collisions. In the threshold model F and a a bs are 
independent and one can fit the data either with small F 
and large a a bs or the vice versa. This can be seen from 
the contour plot of x 2 shown in Figj2j As usual, x 2 is 
defined as, 



X 1 ~ TV 



err(i) 5 



(11) 



where ex(i) and err(i) are PHENIX measurements [13| 
for J/ip suppression as function of collision centrality 
and th(i) is the threshold model prediction. A perfect 
fit is obtained when x 1 < 1. In Fig [2] contours of con- 
stant x 2 /N in F — dabs plane, for fixed values of the 
threshold densities, fm~ 2 and n x —0.36 fm~ 2 

is shown. For a aos varying between 0-2 mb, perfect fit to 
the data could be obtained with F varying between 0.4- 
0.1. Equivalent fit (x 2 /N » 1) could be obtained either 
with F w 0.3, a a bs ~0 mb or with F «0.1 and a aos ~ 
2 mb. The reason is understood also. In the threshold 
model, suppression depends on the transverse density. 
In Fig'121 we have shown the transverse density n p (b, s) 
in Au+Au collisions at impact parameter, b=0,2,4,6,8,10 
and 12 fm. Threshold density for anomalous suppression 
of the state x being small anomalous suppression can 
occur in peripheral collisions also. In contrast, anoma- 
lous J/ip suppression can occur only in sufficiently central 




FIG. 3: Transverse density in Au+Au collisions for various 
values of the impact parameter b=0,2,4,6..(fm) (top to bot- 
tom). The origin is at a distance d = 6/2 from the center of 
the nucleus A. Two solid lines indicate the threshold densities 
for the state J/ip and the state x respectively. 



(b < 6/m) collisions. Unlike for direct J/ip, anomalous 
suppression for the state x occur right from peripheral 
to most central collisions. The suppression is rather like 
nuclear absorption in the Glauber model. In the Glauber 
model, in all centrality ranges of collisions, J/ip's are sup- 
pressed. Naturally, anomalous suppression of the state x 
can approximately compensate for the nuclear absorption 
of J/ip (and x). 

Understanding that both u a bs and F can not be de- 
termined simultaneously from the PHENIX data, for a 
fixed a a bs in the range 0-2 mb, we fit the data by vary- 
ing rij/wj, n x and F. In Fig|2J fit obtained to the data 
for a a b s — 0,1 and 2 mb are shown. Parametric values of 
n J/ipi n x an d F are a -^ so depicted. The threshold model 
with higher state feedback correctly reproduces the cen- 
trality dependence of J/ip suppression in PHENIX mea- 
surements. We note that within the experimental errors, 
J/ip suppression in very peripheral as well as in cen- 
tral collisions is well reproduced even with a a bs=0 mb. 
However, the data require largest contribution from the 
higher state, ^=0.44, when the nuclear absorption is ne- 
glected. The model also reproduces the sudden change in 
the slope above the participant number N part ~ 150. It 
is interesting to note that different parameter sets predict 
nearly identical suppression above N part ~150 . Rather 
the parameter sets are somewhat distinguished in colli- 
sions with Np ar t <150. Apparently, less central collisions 
are more sensitive to the details of the model rather than 
very central collisions. 

To understand better the suppression due to various 
components of the model, in FigJSl we have separately 
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FIG. 4: Solid circles are PHENIX data on the centrality 
dependence of J/ip suppression in Au+Au collisions . The 
lines are threshold model fit to data with parameter sets de- 
picted in the figure. Note that the parameter a J J^ was kept 
fixed during the fit. 



FIG. 5: Contribution of various components of the threshold 
model to J/ip suppression in Au+Au collisions at RHIC are 
shown. The solid circles are PHENIX data [TH on the cen- 
trality dependence of J/ip suppression in Au+Au collisions. 



shown the suppression due to (i) nuclear absorption of 
J/ip (the solid line), (ii) nuclear absorption of J/ip+x 
(the dashed line), (iii) nuclear absorption of J/ip+x and 
the anomalous suppression of \ (the medium dashed 
line) and (iv) nuclear absorption of J/ip+x an d the 
anomalous suppression of J/ip and x (the short dashed 
line). Without any feedback from the state x-> if J/^'s 
are suppressed due to nuclear effects only (<7 a f, s =2 mb), 
PHENIX data in peripheral collisions or in central col- 
lisions are not explained. Suppression is marginally in- 
creased (x-nucleon absorption cross-section is assumed to 
be twice the J/i/>-nucleon absorption cross-section) when 
23% feedback from the state x is included. Suppression is 
further increased when the state x is allowed to undergo 
anomalous suppression (n x =0.27 /to -2 ) in addition to 
the nuclear absorption of x an d J/ip- It can be seen 
that x' s are anomalously suppressed nearly in all the 
centrality ranges of collisions. Anomalous suppression 
of x alone can explain the PHENIX data in peripheral 
and mid-central collisions (N part < 150). Data beyond 
N P art — 150 is underpredicted. Only when J/ip's are also 
suppressed anomalously, PHENIX data throughout the 
centrality range is explained. For nj/^,—3.57 fm~ 2 , di- 
rect J/ip's are anomalously suppressed only in collisions 
with Np ar t > 150. Since, anomalous suppression is re- 
lated to melting of J/ip in deconfined medium, the anal- 
ysis strongly suggests that sufficiently hot (T > 2.1T C ) 
deconfined medium is created in central and mid-central 
{Np ar t > 150) Au+Au collisions. 

Variation of the fitted parameters {nj/^,n x and F) 
with (Tats is interesting. In panel (a), (b) and (c) of Fig[S] 
we have plotted the best fitted threshold densities nj/^, 



n x and the fraction F of the state x- It may be noted 
that the uncertainty in the parametric values are large. 
The quality of the PHENIX data is not good, data points 
are few and the error bars are also large. The uncer- 
tainty in the parametric values will reduce with better 
quality data. It is interesting to note that for a aos =0-2 
mb, the threshold densities for the anomalous suppres- 
sion of direct J/ip and also of the state x is approxi- 
mately a constant. The direct J/ip are anomalously sup- 
pressed beyond a threshold density nj/,/, = 3.57+ 0.17 
/to -2 . Within the uncertainty, the value is consistent 
with rijj^ — 3.73 + 0.29/to -2 , extracted neglecting the 
feedback from the state x [111 • Uncertainty is the thresh- 
old density for the state x is quite large, x states are 
anomalously suppressed beyond n x = 0.32 ± 0.32/m -2 . 
Since in the threshold model threshold density is assumed 
to be proportional to the energy density above which 
the charmonium states get dissociated in QGP, if Tj/^, 

and T x are the dissociation temperatures then ( F ^* ) 4 ~ 

—U± = H.16 ± 11.69. As mentioned earlier, threshold 

n x > 

temperatures of J/ip and x are calculated in the lat- 
tice based models, Tj/^ « 2.12T C and T x w l.2T c [H- 
Within the large uncertainty, the ratio of the threshold 
temperatures extracted from the RHIC data on J/ip sup- 
pression agree with the lattice calculation. 

While the threshold densities Hjm and n x are found 
to be nearly independent of <j a bs, the fraction F is not. F 
linearly decreases with increasing a a bs ■ As discussed ear- 
lier, in the threshold model a a bs a nd F are inter-related 
and one can compensate for the other. Depending on the 
value of dabs, the PHENIX data can be well described 
with F = 0.3 ± 0.1, the maximum (F=0.4) and mini- 
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FIG. 6: In panels (a), (b) and (c) best fitted values for the 
threshold density for direct J/ip, threshold density for the 
state x an d the fraction F of the state \ are shown for various 
values of J /ip -nucleon absorption cross section. The solid line 
in the panel (c) corresponds to the empirical relation Eg |12l 



mum (F=0.2) values corresponding to er a fcs=0 and 2 nib 
respectively. We also observe that an empirical relation, 



F = 0.406 -0.89a abs (fm 2 ), 



(12) 



exist between best fitted F and a a bs- It is also interest- 
ing to note that the maximum value of F «0.4 obtained 
presently ignoring the nuclear absorption agree with the 
fraction F measured in hadron-hadron collisions (25l . [26| , 
where the nuclear absorption is absent by default. 



IV. J/ip SUPPRESSION IN PB+PB 
COLLISIONS AT LHC 

At the Large Hadron Collider (LHC), it is planned to 
collide Pb ions at centre of mass energy y/s—5500GeV. 
Signals of the deconfinement phase transition are ex- 
pected to be better denned at LHC than at RHIC. One 
wonders whether in LHC energy collisions, J/ip will be 
more suppressed than at RHIC energy? The threshold 
model with the parameters fixed from the analysis of 



^™Pb+Pb@LHC (prediction) 
n J/ x I »=(3.74-3.40) fm" 2 
n % =(0.64-0)fnrr 2 ,F=(0.2-0.4) 
°abs=(0-0.23) mb 
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FIG. 7: Threshold model predictions for the centrality depen- 
dence of J/ip suppression in Pb+Pb collisions at LHC energy 
(-y/s^SSOO GeV). The shaded region marks the uncertainty 
in the parametric values of nj/^, n x and F. For comparison 
PHENIX data [HJ on the centrality dependence of J/ip sup- 
pression in Au+Au collisions at RHIC energy are also shown. 



RHIC data can be used to predict the suppression in 
Pb+Pb collisions at LHC. RHIC data on J/ip suppres- 
sion are explained with, nj^ — 3.57 ± 0.17/m -2 and 
n x = 0.32±0.32/m~ 2 , a abs = l±lm& and F = 0.3±0.1 
. nj/^p and n x are independent of energy and their value 
remain unchanged from RHIC to LHC energy. However, 
the parameter a a bs could change from RHIC to LHC en- 
ergy collisions. In SPS energy (^=17 GeV) a a bs ~ 4.4 
mb [H. It is reduced in RHIC energy (^5=200 GeV) col- 
lisions, dabs = 1 ± lmb. At LHC energy (^/s—5500 GeV) 
&abs could be even lower. However, as shown in Figf2]and 
FiglSti in the threshold model, the parameters F and a abs 
are not independent, an empirical relation (Eq fL2)l exist 
between F and a a bs- We assume that at LHC energy 
such an relation exist between a a b s and F and given the 
fraction F, a a b s can be determined using Eq |12l Shaded 
region in FigO show the predicted centrality depen- 
dence of J/ip suppression with nj/^ = 3.57±0.17 fm~ 2 , 
n x = 0.32 ± 0.32 /m 2 and F = 0.3 ± 0.1. At LHC en- 
ergy, for the nucleon- nucleon inelastic cross-section (Inn 
we have used the value o nn—70 mb. In central colli- 
sions, J/ip's are suppressed by a factor of 3-4. For com- 
parison, in Figf7]we have also shown the PHENIX data, 
which as shown earlier are well explained in the thresh- 
old model. Interestingly, threshold model predicts that 
in central Pb+Pb collisions at LHC, J/ip's are similarly 
suppressed (if not slightly less) as in Au+Au collisions 
at RHIC. Reason could be understood as follow: in the 
threshold model, anomalous suppression depends on the 
(local) transverse density. Transverse density depends on 
the nucleon-nucleon inelastic cross- section, and changes 
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FIG. 8: The solid and the dotted line show the transverse 
density in b=0,2,4 and 6 fm Au+Au collisions at RHIC and 
in Pb+Pb collisions at LHC energy. The origin is at a distance 
d = 6/2 from the centre of the nucleus A. 

from RHIC to LHC energy. In FigJSJ the transverse den- 
sity in Au+Au collisions at RHIC and in Pb+Pb colli- 
sions at LHC is compared. In central and mid central 
collisions participant density is less in Pb+Pb collisions 
than in Au+Au collisions. Then for a threshold density 
n J/V" w hile in Au+Au collisions J/tp's could be anoma- 
lously suppressed say in a impact parameter hi collision, 
they may not be anomalously suppressed in Pb+Pb colli- 
sions. Thus J/tp's will be less suppressed in LHC energy. 
Indeed, explicit calculations indicate that under identi- 
cal condition, in central collisions anomalous suppression 
decreases from RHIC to LHC energy. 



V. QGP FORMATION TIME 

QGP formation time is an important parameter in 
ultra-relativistic heavy ion collisions. For example, in 
the Bjorkcn model with longitudinal boost-invariance, 
the energy density is obtained as, 



£bj 



1 dEi 



(13) 



tAt dy 

where r is the formation time, At is the overlap area 
and is the transverse energy per unit rapidity. An 
estimate of energy density produced in a collision then 
depends on the formation time r. QGP formation time 
(r) is unknown. Due to heavy mass, J/tp's are pro- 
duced early in the collisions. Time scale of J/tp formation 
t ~ ~ 0.06 fm, and experimental data on J/tp 

suppression can give information on the QGP formation 
time. 



> 

CD 

O 

— D 

00 
LU 



0.29 npmax 2 
n Pmax 
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FIG. 9: Solid circles are PHENIX data [H on the participant 
number dependence of Bjorken energy density times forma- 
tion time (tEbj) in Au+Au collisions at -^=200 GeV. The 
dashed line is the participant number dependence of the max- 
imum transverse density {n Pmax ) in a Au+Au collision. The 
solid line is the participant number dependence of 0.29rip maa; . 
Empirical relation tebj = 0.29rip max correctly reproduces 
the PHENIX estimate of Bjorken energy density times for- 
mation time. 



PHENIX collaboration have measured the transverse 
energy Et in Au+Au collisions. Since QGP formation 
time is not known, they have tabulated the Bjorken en- 
ergy density times the formation time as a function of 
the participant number. In Fig[9] we have shown the 
PHENIX results for the participant number dependence 
of tebj [33j. tebj increases with the collision centrality. 

In the threshold model local transverse density (see 
EqJSJ) is the key ingredient. If the local transverse den- 
sity exceed the threshold density nj/^p, J/tp's are com- 
pletely suppressed. It is implicitly assumed that the 
transverse density is proportional to energy density and 
in a deconfined matter J/tp's are melted above a thresh- 
old energy density £th- As given in Eq[5l local trans- 
verse density is a function of the impact parameter (b) 
and the transverse position (s). In a collision between 
two identical nuclei at impact parameter b, maximum 
transverse density (n Pmax ) is achieved at the transverse 
position s = b/2. If in a collision J/tp's are anoma- 
lously suppressed, at the minimum n p should exceed 
n J/tpi n Pmax ^ n J/tp- I n Fig[5]the dashed line show the 
participant number dependence of the maximum trans- 
verse density n„ = n„(b.s = b/2). n„ increases 
with collision centrality. The solid line in Fig[5] indicate 
that empirically, in all the centrality ranges of collisions 
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TABLE I: Formation time r for a choice of critical tempera- 
tures for two and three flavored QGP. 





flcnrrvrorl HPP 


O-lldVOI fc!U V^Vjr_r 




9qgp—o( 


n r ^ ^ ^ — AT t\ 


i c (ivie v ) 


T l rm J 




i ao nn 

14U.UU 


U.ol It U.Uo 


n O/i -i- n no 

U.Z4 It U.UZ 


i Kf\ nn 
loU.UU 


n O/i -i- n no 
U.Z4 it u.UZ 


n io i n no 
U.lo zh U.Uz 


i fin nn 


n ic-Lfi no 

U.lo zt U.UZ 


n 1/1 -i- n ni 

U. 14 zt U.Ul 


i vn nn 
1 / u.uu 


n i a -\- n ni 
U. 14 zt U.Ul 


nil i n ni 
U.ll zt U.Ul 


i sn nn 

loU.UU 


nil _i_ n m 
U. 1 1 zt U.Ul 


n no 4- n ni 
u.uy zt u.ui 


190.00 


0.09 ±0.01 


0.07 ±0.01 


200.00 


0.07 ±0.01 


0.06 ±0.01 


210.00 


0.06 ±0.01 


0.05 ±0.00 


220.00 


0.05 ± 0.00 


0.04 ± 0.00 


230.00 


0.04 ±0.00 


0.03 ±0.00 


240.00 


0.04 ±0.00 


0.03 ±0.00 



0.09(0.07) fm for 2(3) flavored QGP. Formation time r 
is rather small compared to initial time Ti=0.6 fm ob- 
tained from hydrodynamic analysis of experimental data 
in Au±Au collisions at RHIC [35]. However, we note that 
the formation time scale 0.09(0.07) fm is rather consis- 
tent with Color Glass Condensate (CGC) models 36 . 

Before we summarise our results, we would like to 
mention a deficiency of the model. Transverse expan- 
sion of the system is neglected, it is assumed that all 
the J ftp's are suppressed before the transverse expansion 
sets in. However, hydrodynamics model simulations of 
Au±Au collisions at RHIC indicate that even the QGP 
phase undergoes certain amount of transverse expansion. 
With transverse expansion, transverse density is diluted 
and J ftp's which are anomalously suppressed in a static 
medium may not be suppressed in an expanding medium. 
Thus the threshold density nj/^ = 3.57± 0.17/ra -2 rep- 
resent an upper limit to threshold density. The QGP for- 
mation time as tabulated in tableU could be even lower. 



PHENIX data 33] on tebj exhibit a simple relationship 
with the maximum transverse density n Pmax , 



VI. SUMMARY AND CONCLUSIONS 



tebj = 0.29n 



Pmax 



(14) 



It is an important observation. We have checked that 
the relation t£r r oc ?i„ is consistent with measure- 

Pmax _ 

ments of tebj in other energies e.g. \fs— 130 GeV and 
19.6 GeV also. Since n p > nj/^, empirical relation 
(EafU)) imply, tebj < 0.29n 2 ^. If the threshold energy 
density for J ftp dissociation is Eth then it is necessary 
that, 



0.29n 



> TE B J > TETH 



(15) 



In otherwords, the threshold density in the threshold 
model give the upper limit of the threshold energy den- 
sity for J ftp dissociation times the QGP formation time. 
Assuming that the J ftp's are melted above a temperature 
Tth = 2.1T C , given a critical temperature, the standard 
relation between energy density and temperature, 



9QGpt^T th 



g QGP — (2.1T c )\ (16) 



can be used to compute the upper limit of the QGP for- 
mation time. 



In table. HI for threshold density 



l J/tp 



= 3.57 ± 0.17 



/m -2 , we have tabulated the upper limit of the formation 
time t for a two and three flavored QGP. QGP formation 
time is low. For critical temperature T c — 140-200 MeV, 
it varies from 0.07(0.06)-0.31(0.24) fm. Critical temper- 
ature for the connnement-deconfmcmcnt transition has 
been accurately estimated in recent lattice calculations, 
T c = 192±7±4 MeV As seen from tableffl T c - 192 
MeV corresponds to formation time r ranging between 



To summarise, in the QGP motivated threshold model, 
we have analyzed the mid-rapidity PHENIX data on the 
centrality dependence of J ftp suppression in Au±Au col- 
lisions. In the threshold model, in addition to the normal 
nuclear absorption, J/tp's are anomalously suppressed, 
such that, if the local transverse density exceeds a thresh- 
old value rij/^, all the J/tp's are absorbed. We have ex- 
tended the threshold model to include the feedback from 
higher states x an d The extended model has four 
parameters, the threshold density {nji^f) for anomalous 
suppression of direct J ftp, the threshold density (n x ) for 
anomalous suppression of the state %, the fraction F of 
the states \ and the J/i/)-nucleon absorption cross section 

(fafcs)- 

We observed that the centrality dependence of J ftp 
suppression can not determine the fraction F of the state 
X and the nuclear absorption cross-section a a bs indepen- 
dently. In the threshold model they are inter-dependent. 
Anomalous suppression of the state x, to certain ex- 
tent, can mimic the nuclear absorption of J ftp and x 
and the PHENIX data can be fitted either with large F 
and small a a bs or the vice- versa. For a aos =0-2 mb, we 
have fitted the PHENIX data by varying Ujm, n x and 
F. Data are well described for nj/^ = 3.57 ± 0.17 /m 2 , 
n x — 0.32 ± 0.32 fm 2 . The uncertainty in the paramet- 
ric values could be reduced if the quality of the data is 
improved. The fraction F depend on the value of a a b Sl 
F = 0.3 ±0.1, the upper and lower values corresponds to 
o"a6s = and 2 mb respectively. 

Threshold model with parameters fixed from the 
Au±Au collisions is then used predict for the central- 
ity dependence of J ftp suppression in Pb±Pb collisions 
at LHC energy. At LHC energy the model predicts that 
in central collisions J/tp's are suppressed by a factor of 
3-4. Suppression pattern is similar to that in Au±Au 
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collisions at RHIC. In central collisions, possibly J/ip's 
are slightly less suppressed in Pb+Pb collisions at LHC 
energy than in Au+Au collisions at RHIC. 

We have also estimated the QGP formation time. On 
the basis of the empirical observation that in all the cen- 
trality ranges of Au+Au collisions, the square of the max- 
imum transverse density times a factor 0.29, correctly re- 
produces the PHENIX measurements on Bjorken energy 
density times the formation time, we have argued that 
the threshold density gives the upper bound for the QGP 



formation time. For critical temperature T c m= 192 
MeV, QGP formation time is estimated as t w 0.09(0.07) 
fm. 

To conclude, threshold model analysis of J ftp suppres- 
sion in Au+Au collisions at RHIC, strongly suggest that 
in central and mid-central collisions (N par t > 150) QGP 
is produced at sufficiently high temperature T > 2.1T C , 
to melt the directly produced J ftp's. The analysis also 
indicates that QGP is produced rather early in the colli- 
sions. 
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